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Abstract. We study the flashing ratchet model of a Brownian 
motor, which consists in cyclical switching between the Fokker- 
Planck equation with an asymmetric ratchet-like potential and 
the pure diffusion equation. We show that the motor really per- 
forms unidirectional transport of mass, for proper parameters of 
the model, by analyzing the attractor of the problem and the sta- 
tionary vector of a related Markov chain. 



1. Introduction 

Nano-scale or molecular devices which use energy but not momen- 
tum to generate transport are called Brownian motors. Such phenom- 
ena arise in various areas of science, from intracellular transport to 
nanotechnology [2j El EJ El E] . 

The general relation for various types of fluctuation- driven motors 
looks like [2] 

(1) Pt = o-p xx + {^ x p) x , x <E (0,l);t> 0. 

Here p is the unknown density, a is the diffusion coefficient, and 
ty(x,t) is the potential. For the flashing ratchet, an autonomous po- 
tential ip is switched on and off cyclically [2], i.e. [71 H] 

(2) *(x,t) = h(t)i/>(x), 
where 



(3) h(t) 



1, nT < t < nT + T tr , n = 0, 1, . . . , 
0, nT + T tr < t < nT + T, n = 0, 1, 
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A typical ratchet-like potential ip with k teeth, k > 1, is l//c-periodic 
in x and has a unique local (and, hence, global) minimum within each 
period. 

In [7J it was shown that the behaviour of the flashing ratchet system 
(with Neumann boundary conditions) is in some sense close to the 
behaviour of a certain Markov chain. It was observed that having 
this at hand it is possible to verify transport via comparing eventual 
distribution of mass between the "wells" of the potential ip, i.e. the 
line segments with end points at successive maxima of ip(x). Any 
inequality in this distribution would mean transport. In particular, it 
was shown that exactly this takes place for proper parameters of the 
ratchet and k = 2. However, the proof was not completely consistent, 
being based on incorrect time asymptotics of the second derivative of 
the Green function (a power function instead of an exponential one). 
Generalization of this claim to the case k > 2 was mentioned as an 
open problem in jH [7] . 

In this paper, we give evidence of unidirectional transport for any k > 
1. Namely, we show (Theorem I2.2p that for certain parameters of the 
flashing ratchet, after a sufficiently large number of cycles the amount 
of mass in the wells of the potential is strongly decreasing/increasing 
from the left to the right, provided the minima of ip(x) are located in 
the left/right halves of the wells. 

It is important to note that the transport provided by the flashing 
ratchet is due to flashing fl3]) only, since both pure diffusion {h = 0) and 
"perpetual ratchet" {h = 1) with a periodic potential rapidly approach 
their equilibria without any specific right or left drift tendency. 

Let us also recall that there is a connection (see e.g. [TJ EJ [7J) be- 
tween the flashing ratchet, especially the fact that it produces unidi- 
rectional transport, and Parrondo's paradox in game theory, where a 
well-scheduled alternation in playing two fair or losing games becomes 
a winning strategy. 

The paper is organized as follows. In the next section, we present 
the problem more rigorously, give necessary notations and facts, and 
formulate the main result (Theorem I2.2p . In the third section, we 
demonstrate that the so-called discrete ratchet, which generates the 
Markov chain, behaves in a way similar to the claimed behaviour of 
the flashing ratchet. The proof of the main result is provided in the 
final section. 
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2. Preliminaries 

We consider the boundary value problem for the flashing ratchet 
equation with Neumann boundary conditions: 



(4) 



Pt = <?pxx + h(t)(bp) x , x e (0, 1); t > 0, 

ap x + h(t)bp = 0, x — 0, l;t>0, 

i 

p > 0, / p(x,t)dx = 1, t > 0. 
o 



Here 6 is the x-derivative of the potential ip and h is given by 
The ratchet phase time periods are of length T tr , and the pure diffusion 
periods are of length Tdiff = T — T tr . We denote 

T = vTdiff- 

Following [5], the potential ip(x) is assumed to be a C 4 -smooth 
function on [0, 1] of period 1/k, with k > 1 being a fixed integer, having 
maxima at points and minima at points a, and being monotone 
between these points (a ratchet-like form), where 

i — 1 

(5) Xi = — — , z = 1, . . . , k + 1, 

k 

(6) cii = a + Xi, i = 1, . . . , k. 

The positive parameter a should be less than 1/k. 

The problem can be completed with the initial condition 

(7) p(x,0) =p Q (x), xe{0,l), 
such that 

i 

(8) p (x) > 0, / Po (x)dx = 1. 





The existence of a periodic orbit for is provided by 
Theorem 2.1. (see [3 Theorem l]j Assume that 

(9) 2tt 2 t - AT tr > In 2, 

where X is a certain constant depending only on the potential (see [7] 
for its exact value). Then problem fll]) has a unique T -periodic solution 
P s - 

It is also known [71 H] that the periodic orbit p s eventually attracts 
all the solutions p of (HJ), namely, 

(10) lim \\ P (;t n )-p S (;t n )\\ H2 =0. 

n— >oo, t n =t+nl 
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For brevity, sometimes we will write simply p s or p s (x) for the function 
p s (x,0) = p s (x,nT). In [7J it is shown that 

(11) \\f\\ H ><Ro = ^L±±. 
Denote 

Xi+i 

pi = J p s (x)dx, i = l,...,k. 

Xi 

The main result of the paper is 
Theorem 2.2. If a < h-, then there exist a, T tr and T such that 

(12) p\> &>■■■> p s k . 

Theorem 12.21 means, in particular, that, given any initial distribution 
of density, after a large number of cycles there will be more mass on 
the left than on the right. 

Remark 2.1. If a > -h-, then Theorem 12.21 implies p\ < p s 2 < ■ ■ ■ < p s k , 
to see this it suffices to make the change of variables x — > 1 — x. 

The discrete ratchet acts a follows. During the ratchet phase it 
simply concentrates all the matter from any segment [x^Xj+i] at the 
point aj. Thus, if 

p* = J po(x) dx, 

X i 

then at the moment T tr the density becomes 

k 
i=l 

During the diffusion phase we have the same diffusion as for the flashing 
ratchet. Then this process is repeated periodically. 

Denote by d the Wasserstein metric (see e.g. [B]) on the space of 
probability measures on [0, 1]. For a continuous function / and a prob- 
ability measure p on [0, 1], we use the bra-ket notation as follows: 

i 

(P> /) = / / d P- 
o 

The convergence in Wasserstein metric implies *-weak convergence of 
probability measures, i.e. 

(13) d(p n , p) -> (p n - p, f) -> 0, / G C[0, 1]. 
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At the end of the ratchet phase, it is possible to estimate the distance 
between the solution to flU),© and the outcome of the discrete ratchet: 

Lemma 2.1. (see [7J Corollary 3}) Let p be a solution to flU),© with 
\\Po\\h 2 < Ro- Then for sufficiently large T tr (the lower bound on T tr 
depends on the potential only) one has 

( k \ 2 1 2 T 

(14) d\p(-,T tr ),J2^5 ai j < R^(l +Cl )^^+mm{C,ae XT ^ 2 ,l}. 

The constants c\ and C\ depend on the potential only. 

Let us describe how the discrete ratchet generates a Markov chain. 
Consider the heat equation with Neumann boundary conditions: 



Vs = Vxx, x e (0, 1); s > 0, 
y x = 0, x = 0, l;s>0. 



= ex P(-g/ 4g ), 

2\/TXS 



(15) 
Let 



and 

oo 

(16) G{x,s) = T s (x + 2n). 

n=— oo 

Note that G is 2-periodic and even in x, and G(l + x,s) = G(l — x,s). 
The Green function for (fl5|) is [TO] 



n _s g(£, x, s) = 1 + 2 cos(n7r^) cos(n7ra;) exp(— n 2 7r 2 s) 

(1'J n=l 

= G , (x + e,s) + G(x-^s). 
Now introduce the following matrix: 

Xj+l 

(18) P = P(t) = (pij), pij = J g(ai,x,r)dx. 

Since the initial distribution of mass between the segments [xi, Xi+±] 
is given by the vector /x* = (//*), the outcome of the discrete ratchet at 
the moment T will have the distribution of mass between the segments 
described by the vector /x*P(r) (cf. [7]), at 2T it will be /i*P 2 , and so 
on. 
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3. The stationary vector of the Markov chain 

In order to prove the main theorem we need first to study the sta- 
tionary vector of the Markov chain generated by the matrix P. 

We recall that an m x m-matrix with positive entries is called ergodic 
if the sum of the elements in every row is equal to one. The eigenvalue 
1 of any ergodic matrix A is simple, and there exists a unique vector £ 
satisfying 

m 

(19) &>0, i = l,..., m , ]T& = U = £P. 

i=i 

We will call it the stationary vector of A. Let us also introduce the 
number 

k(A) = min \yA — y\. 

m 

E S/i=0,|2/|=l 
1=1 

Note that, for an ergodic matrix A, this number is positive, and k(A) — >■ 
1 as all the elements of A approach 1/m. 

The matrix P is ergodic. Denote by fi s = (/i|) its stationary vector, 
which is also the stationary vector of the corresponding Markov chain. 
The following result holds. 

Theorem 3.1. For a < ^ and r large enough, there is a constant 
c > independent of r such that 

(20) £ > i4 + ce~« 2 \ a4 > /4 + ce'^, . . . , ^ > ^ + ce'^. 
Its proof requires 

Lemma 3.1. Let A = (a^) be an ergodic m x m-matrix satisfying the 
following conditions: a) for any column (say, j-th, j < m) there exists 
a number s = s(j) so that one has < a«j+i provided i > s, and 
a ij > a i,j+i provided i < s, b) there exists a constant d > such that 

A 1 >A 2 + d, A 2 >A 3 + d,..., A m _ x >A m + d, 

where Aj is the sum of elements in the j-th column. Then the stationary 
vector £ of A satisfies 

(21) 6 > 6 + Md, 6 > 6 + Md, . . >U + Md, 

where M is the minimum of the elements of the last (i.e. m-th) column. 
Proof. Consider the set 



B=lye 



^Vi = l,yi>y2 + Md, 2/2 > 2/3 + Md, 



i=l 
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. 2/m-l > Vm + Md > Vm > M 



This set is compact and convex. Moreover, B is invariant for the map 
A : y (->■ yA. In fact, let y E B. Then 



m m 



Y Y Y 11 • Y "i* = '• 

i=l j=l j=l i=l 

Fix any I = 1, . . . ,m — 1. Then 

m m 

- = Yy^ au ~ a ^+0 - ^^(q( q ^ - a M+o 

i=l t=l 

= y s {i){Ai - A l+1 ) > y s{ i)d > Afd. 

Finally, 

i=l i=l 

By Brouwer's fixed point theorem, ^4. has a fixed point in B, which 
should coincide with the stationary vector. □ 

Proof. (Theorem 13. ip Let us show that P satisfies the conditions of 
Lemma I37T1 Set s = [^p] (the integer part) for any j. Firstly, let us 
check if pij > Pij+i for i < s. Consider the function 

4>(y) = / g(ai,x,T)dx, <y < 



k 
y 

It suffices to show that it is decreasing. Observe that a* < | and 
cos(7raj) > since i < s. Now, 

<j>'(y) = g(ai, y + l/k,r)- g(a h y, r) 

oo 

= 2 cos(727raj) (cos(mn/ + rrn/k) — cos(nny)) exp(— n 2 7r 2 r) 

n=l 

= 2 cos(7rcij)(cos(7n/ + ir/k) — cos(7n/)) exp(— 7r 2 r) + o(e _7r2T ) 

< 2cos(7rai)(cos(7r/A;) - 1) exp(-7T 2 r) + o{e~^ T ) < 
for large r. 

The claim that p^ < Pij+i for i > s can be proven similarly taking 
into account that a« > \ and cos(7raj) < for i > s. 
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Without loss of generality (i.e. for large r) we may assume that the 
function G(x, r) is decreasing in x on the segment [a, 1 — a], and 

(22) G x (x,t) < -Ce~ n2T ,a <x < 1 -a, 

with some constant C > 0. Really, from ( IT71) we get the following 
representation: 

j OO 

(23) G(x, t) = - + cos(n7rx) exp(— n 2 n 2 T). 

n=l 

Thus, 

oo 

G x (XjT) = — Tj n7r sin(r27ra:) exp(— n 2 ix 2 r) 

n=l 

= — 7T sin(7rx) exp(— 7T 2 r) + o(e _7r2T ) 

< -7rsin(7ra) exp(-7r 2 r) + o(e^ T ) < -Ce _ ^ T . 
Take any / = 1, . . . , m — 1. We have to see that 

(24) > P l+1 + d 

with some d independent of I. Consider the function 



k v ^ 

¥{y) = ^2 J 9{<H,x,t) dx, < y < 



i=i 
y 

Then 



fc-1 



9? 



'(y) = Y)g(a + l)/k, y + 1/k, r) - g(a + (i - l)/k, y, r)} 



k 



^2[G(a + y + i/k, t) + G(a - y + {i - 2)/k, r) 



i=l 



-G(a + y + (i - l)/k, r)-G(a-y+(i- l)/k, r)] 

= G(a + y + l,T) + G(a-y-l/k,T)-G(a+y,T)-G(a-y + (k-l)/k,r) 

= G(l-a-y,r)+G(y-a+l/k,r)-G(a+y,r)-G(a-y+(k-l)/k,r). 

The length of the segments [a + y, y — a + 1/k] and [a — y + (k — 
l)/k, 1 — a — y] is r — 2a. Therefore (122"]) implies 

(25) ^'(y) < C(4a - 2/k)e-^ T . 
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Thus, (1Mb holds with 



^_ C(-4a + 2/k)e-^ T 
d ~ k ' 

It remains to observe that due to (fTTD one has 

(26) 



M = \ + 0{e~^ T ) 



for the minimum of the last column of P. 



□ 



4. Proof of the main theorem 
Proof. (Theorem 12.21) Observe that there exist sequences 

T tryU -)■ oo, r n -)■ oo, o n ->■ 

satisfying 

(27) o- n e AT ^ n/2 ->■ 
and 

(28) 27r 2 r n - AT tr>n > In 2. 

Let T n = T tr n + — . Then we can find the corresponding T n -periodic 
solutions p s,n , and introduce obvious notations p s ' n , p s ' n etc. for the 
corresponding values. It suffices to show that 

(29) \p s ' n -p s ' n \ = o{e-^ Tn ) 

as n — > oo. In this case, setting r = r n ,T tr = T tr>n and T = T n for n 
large enough, we would get (fl2"|) from (I20I) . 



Let 



p* s ' n (x 



t) = } J Pi' n g(ai, x, a(t - T 4r ,„)), t > T tr , n - 



i=l 



Then 



\p s r - {p s ' n P{r n )) 3 



Xj + 1 



J p s > n (x } T n )-p* s > n (x,T n )dx 
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(30) 



x j + i I k 

I (p s ' n (-,T tr )-Ep-'X,i 

Xi \ i=l 



+ 2 ( p s ' n (-,T tr ) - EPi' n ^,cos(7r-) )cos(7ra;)e 



i=i 



+ 2 J2 (p s,n (-, T tr) ~ EPi' n( ^,cos(m7r-) ) cos(mvrx)e- m 77 T " dx 

m=2 \ i=l 



Due to (TT3D and (H 



j = 1, . . . ,n. 

The first summand is zero, the third is o(e 
with p(x,t) = p s ' n (x,t), p* = p s ' n , one has 

P S '"(-,^)-gA S 'X,cos(7r-)^0, 

so the second summand from ( 1301) is also o(e _7r2r ™) as n — > oo. Thus, 
(31) 



\p s ' n - p s ' n P(r n )\ = o(e~* Tn ). 



It remains to observe that 



(32) 

and n(P(T n )) -> 1. 



|p s ' n - p s ' n P(r n 
«(P(t-„)) 



□ 
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